Hard thermal effective action in QCD through the thermal operator 
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We derive in a simple way the well known hard thermal effective action for QCD through the 
application of the thermal operator to the zero temperature retarded Green's functions. This deriva- 
tion also clarifies the origin of important properties of the hard thermal effective action, such as the 
manifest Lorentz and gauge invariance of its integrand, by relating them directly to the properties 
of the corresponding zero temperature effective action in the hard regime. 
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I. INTRODUCTION 

The high temperature limit of thermal QCD is of much 
interest because of its relevance in the study of the phys- 
ical properties of the quark-gluon plasma [J 0] • Several 
important features in this regime are embodied in the 
leading hard thermal effective action, whose overall coef- 
ficient is proportional to T 2 , where T denotes the equi- 
librium temperature. Such leading contributions to the 
effective action arise at one loop from diagrams where 
the internal momentum is hard (being of the order T) 
and is much larger than any external momenta. The 
hard thermal effective action in QCD, which enjoys var- 
ious symmetry properties (at the integrand level) such 
as manifest Lorentz and gauge invariance OjLljk has also 
been studied from various points of view [ETEol IS & El • 

The purpose of this note is to present a simple deriva- 
tion of this action which may explain the origin of these 
interesting symmetry properties. The derivation is based 
on an interesting relation between Feynman graphs at 
finite temperature and the corresponding zero tempera- 
ture graphs [111) 03 > which holds both in the imaginary 
time as well as in the real time formalisms [jj, [2j, LL3( • This 
relation, known as the thermal operator representation, 
arises as a consequence of the fact that the thermal prop- 
agator for a bosonic field can be related to the zero tem- 
perature one through a simple thermal operator which 
carries the entire temperature dependence and has the 
explicit form 



0^(E) = l + n B {E)(l-S{E)). 



(1) 



Here E = + m 2 ,S(E) is a reflection operator that 

takes E — > — E and Ub{E) represents the bosonic distri- 
bution function. (For a fermionic field a similar relation 
holds with tir — > — rip [12].) This relation between the 
finite temperature Feynman graphs and the zero tem- 
perature ones is calculationally quite useful and allows 
us to study directly many questions of interest such as 
the cutting rules at finite temperature [Til ] . Furthermore, 
the relation between the retarded thermal Green's func- 
tions and the forward scattering amplitudes for on-shell 
thermal particles [HI has been clarified through the ap- 
plication of the thermal operator representation to the 



corresponding zero temperature forward scattering am- 
plitudes [l6j . In this paper, we further demonstrate the 
simplicity and the utility of the thermal operator repre- 
sentation by deriving the hard thermal effective action for 
QCD from the zero temperature retarded amplitudes in 
the hard regime. This derivation directly associates the 
origin of the interesting properties of the hard thermal 
effective action, such as Lorentz and gauge invariance of 
the integrand, to those of the zero temperature effective 
action. 

The paper is organized as follows. In section II, we 
show that in the hard region at zero temperature, the 
forward scattering amplitudes are Lorentz and gauge co- 
variant and obey simple Ward identities. These features, 
together with the fact that the leading one loop contri- 
butions are quadratic in the hard internal momentum, 
are sufficient to determine uniquely all the hard n-point 
gluon amplitudes in terms of the hard gluon self-energy. 
An explicit example of how this works for the 3-point 
gluon amplitude is discussed in more detail in the ap- 
pendix. In section III, we are thus able to write down 
a generating functional (effective action) for all the hard 
n-point gluon amplitudes at zero temperature. Through 
the application of the thermal operator, it is then imme- 
diate to arrive at the well known form of the hard thermal 
effective action. In this way, the symmetry properties of 
the action, such as Lorentz and gauge invariance, can be 
directly understood in terms of the properties of the zero 
temperature retarded amplitudes. We conclude with a 
brief summary in section IV. 



II. HARD FORWARD SCATTERING 
AMPLITUDES AT ZERO TEMPERATURE 



The forward scattering amplitude associated with the 
retarded gluon self-energy at one loop, II™, (p), can be de- 
scribed as in Fig. 1, where "R" denotes a retarded prop- 
agator while the cut line with momentum k represents 
an on-shell particle scattering in the forward direction 
[l^ . In the hard region where the internal momentum 
k ^ p, with k 2 = 0, we can expand the denominator (of 
the retarded propagator where ko is to be understood as 



2 



k 



in the graph as 



(k+p) 2 2k -p (2k -pY 



(2) 



The important thing to note here, for later use, is that 
this expansion does not involve any factor of p 2 in the 
denominator. 
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FIG. 1: Forward scattering graphs for the retarded gluon self-energy at one loop. Wavy and dashed lines denote respectively 
gluon and ghost particles. Color and Lorentz indices are suppressed for simplicity. 



Similarly, expanding the numerator of the graph in pow- 
ers of the external momentum, we find that the lead- 
ing hard contribution to the gluon self-energy, I¥^ v (p) = 



ten as I If 



in the SU(N) Yang-Mills theory can be writ- 



IV Op) 



2g 2 N 



d^k 
(2^)3 



S+(k 2 



kuPv ~f k u pn 2 kuk v 

+ p 



k ■ p 



(k-pf 



(3) 



where we have defined 



6+(k 2 ) = 6(k )6(k 2 ). 



(4) 



There are several things to note from the structure of 
the gluon self-energy in (|3|). The quantity in the paren- 
thesis is manifestly Lorentz covariant being of zero degree 



in the internal momentum. Consequently, the integrand 
is quadratic in k for large values of the internal momen- 
tum and, as we will show shortly, this leads to the leading 
contribution of order T 2 at high temperature. To apply 
the thermal operator, we need to carry out the integra- 
tion over ko [12| in which case (J3J) takes the form 



IV (p) 



2g 2 N 



d 3 k 1 
(2tt)3 2E k 



kjxPv + k u p^ 
k ■ p 



2 k^ky 
9 (k-pf 



(5) 



r 



where £7. = \k\ and we have defined k^ — (1, — k). (We 
use a metric with signature (+,—,—,—).) ft is worth 
remarking here that the self-energy in ([5]) is gauge inde- 
pendent and satisfies the Ward identity 



^iV(p) = 0. 



(6) 



Next, let us consider the retarded n-point gluon for- 
ward scattering amplitude at one loop. As discussed in 
[lfl, this is described by graphs with a single on-shell 
cut propagator, all other propagators in a graph being 
retarded/advanced. Thus, one can write such an ampli- 
tude in the form (compare with P|)) 



T~iaia2---aTi 



d 4 fc 
(2^)3 



6 + (k 2 



(7) 
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where 7 represents the forward scattering amplitude for 
an on-shcll particle of momentum k. For large values 
of k, we can again expand the denominators in 7 (see 
@) in powers of p 2 /(2k ■ p), where p denotes a linear 
combination of the external momenta. Furthermore, we 
can also expand the numerator in powers of pi^/\k\. The 
first term in the expansion comes from the leading term in 
each of the retarded propagators of the form l/(k-p) and 
a numerator independent of pi. The contributions to 7 
from these superleading terms, which are of degree one in 
k, however vanish by symmetry under k — > — k. The next 
term in the expansion which does provide the leading 
quadratic contribution to the n-point gluon amplitude 
comes from terms in 7 which are functions of degree zero 
in k [H. 

In the hard regime, to leading order, these amplitudes 
obey simple linear Ward identities. This follows as a con- 
sequence of the fact that diagrams with external (open) 
ghost lines, which appear in the BRST identities, have 
one less power of k in the numerator compared with 



the corresponding diagrams involving only external gluon 
lines. As a result, the first term in the expansion of such a 
graphs is also naively quadratic in k much like the lead- 
ing contribution for the gluon amplitude. However, in 
the case of graphs with external ghost lines, these lead- 
ing contributions cancel by an eikonal identity when all 
graphs are added. Therefore, the contribution of the am- 
plitude becomes subleading compared to the gluon am- 
plitude. In the appendix, we illustrate how such a can- 
cellation takes place for the one loop ghost-ghost-gluon 
(3-point) amplitude. As a consequence, the leading terms 
in the n-point gluon amplitude obey simple Ward iden- 
tities. For example, to leading order, the 3-point gluon 
amplitude satisfies the relation 

P^f u c x ( Pl ,p 2 , P3 ) = igf abc 0V( Pl ) - IV(p 2 )] , (8) 

showing that it is related to the hard self-energy. Simi- 
larly, the Ward identities relate the 4-point gluon ampli- 
tude to the 3-point amplitude as 



J 



P f Ft v %{PUP2,PZ,PA) = ig [f Cde Tf v \{PUP2,PZ+PA) + f Me r™\(puP2+p i ,P3) + f^^fa + P 4,P2,P 3 )] ■ (9) 



The above properties are sufficient to determine uniquely 
to leading order, all the n-point gluon amplitudes in 
terms of the hard gluon self-energy at zero temperature. 
This is discussed in detail in the appendix for the 3-point 
gluon amplitude, but similar arguments hold for all n- 
point gluon amplitudes. (We would like to emphasize 
here that, in general, the solution to the Ward identi- 
ties such as (|8]) is not unique since the transverse part 
of the amplitude is not fully determined by the identity. 
However, in the hard region, the fact that the leading 
term has an integrand of zero degree in k and that the 
denominators have the form k ■ p as in ([2]) is sufficient to 
determine all the amplitudes uniquely.) 



III. THE HARD EFFECTIVE ACTION AND ITS 
PROPERTIES 

The simple Ward identities satisfied by the hard n- 
point gluon amplitudes can be written in a compact form 
as 

%b 5T[A] , . , abcAC , 5T[A] 



D al 

SAUx) 



= (S ab d^ - gf abc A 



= 0, (10) 



where T[A] denotes the generating functional (effective 



action) and D?f is the covariant derivative. Relation (fTO)) 
implies that the effective action, T[A], is invariant under 
an infinitesimal non-Abelian gauge transformation with 
parameter uj a (x), namely, under 



A a 



Aa(uj) 



Dfuj b , 



T[A] -» r[A(" 



(11) 



the infinitesimal change in the effective action is given by 



6T[AM] 



Suj (1 



5A 



6(u0 



5T[A] 



j=0 



Scu a 5A b „ 



- d ? 6 -JaJ = °> (12) 



where we have used a compact notation suppressing all 
the intermediate integrations. 

Using J7]) and performing the kg integration, we can 
write the effective action to leading order in the form 



r[A] 



d 4 x 



d 3 k 1 
(2irf 2E k 



l(A,k). 



(13) 



From our earlier discussion, we note that 7 has the fol- 
lowing characteristics. 

1. It is gauge invariant. 

2. It is a Lorentz invariant function of zero degree in 
kfj, with ka = \k\. 

3. The integrand involves denominators with products 
of factors of the form k ■ p where p is some linear 
combination of external momenta. 

As we have discussed earlier, in the hard region, all the n- 
point gluon amplitudes are determined uniquely from the 
hard gluon self-energy. Therefore, any effective action 
satisfying the properties listed above and yielding the 
correct one loop hard gluon self-energy would correspond 
to the unique hard effective action. Comparing with the 
hard gluon self-energy in ©, it follows that the hard 
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effective action at zero temperature can be written in 
the form 



T[A] 



0?J d4x jSr k ^Mm^M), (14) 



where the gauge covariant potential, introduced in [icl |. 
has the form 



A a Jx,k) 



k-D 



-k v F 



A" 



" ^k-D^ 



d^k-A) , (15) 



and F^ v denotes the non-Abelian field strength tensor 

F; v = d^Al - d v Al + gf abc AlAt. (16) 

The gauge covariant non-Abelian potential, A^x, k), is 
in general a nonlocal function, where the nonlocality 
arises from retarded line integrals along the direction of 
k. 

Once the hard effective action (fT4")l is determined at 
zero temperature, the hard thermal effective action can 
be determined by applying the thermal operator in the 
following manner, ft is worth recalling [12, UH that the 
thermal operator acts on functions of energy in the inte- 
grand before the spatial momentum integrations are car- 
ried out. Furthermore, in spite of the fact that the per- 
turbative expansion of the covariant potential in (|15[) in- 
volves retarded propagators, the thermal operator leaves 
such propagators invariant 16] so that the covariant 
gauge potential is unaffected by the application of the 
thermal operator. As a consequence, the only term in 
the integrand on which the thermal operator {l} acts 
nontrivially is the energy denominator (coming from the 
on-shell cut propagator) 

° {T){Ek) 2^E k = YE k (1 + 2riB ^')) ■ 

Thus, the application of the thermal operator immedi- 
ately leads to the temperature dependent hard thermal 
effective action as 

T^)[A] = 0^j d± x J d*k?^A°(x,k)A^(x,k). 

I ' (18) 

Since E k = |fc|, the radial momentum integration can be 
carried out using the standard integral (with the Boltz- 
mann constant set to unity) 



dk kriB(Ek) = I dk 



tt 2 T 2 



e k ' T -l 6 
which yields the hard thermal effective action 



(19) 



T m [A] = 9_2_K J d , x J ^ A a (x>k)A , a{x ^ (20) 



Here d£l represents the angular integration over the unit 

vector k. Finally, integrating this expression by parts, we 
can recast it in the well known form of the hard thermal 
effective action 0, H| 



r W[ A ] 



d 4 x 



dn 

47T 



/ 



a b 



(k-D) 



(21) 

= g 2 T 2 N/Q as the square 



where we have identified 
of the thermal gluon mass. 

It is now straightforward to extend the above result 
for the pure Yang-Mills theory to QCD with quarks in 
the fundamental representation. To this end, we recall 
that the thermal operator representation works also for 
theories involving fermions [l 21 ] . In this case, as we have 
pointed out earlier, the thermal operator, relating the 
propagator at zero temperature to that at finite temper- 
ature, has a form similar to the one given in ©, with 
ub{E) — * —np(E), where tif denotes the fermionic dis- 
tribution function. Because of the linearity of the Ward 
identities in the hard regime, the contribution of the hard 
quark loops to an amplitude, which is additive, indepen- 
dently satisfies simple Ward identities like those given in 
Eqs. ([5]) and ©3- As shown in the appendix, these 
identities together with the fact that the leading order 
terms in the integrand are Lorentz covariant functions of 
degree zero in the internal loop momentum arc sufficient 
to determine uniquely all the higher n-point gluon ampli- 
tudes in the hard regime in terms of the corresponding 
gluon self-energy. Given the above properties of hard 
amplitudes as well as the transversality condition in ^ , 
it follows that the integrand for the self-energy with a 
hard quark loop at zero temperature necessarily has a 
structure similar to the one given in ([3]), upto an overall 
multiplicative factor which can be calculated easily. Con- 
sequently, adding this contribution to that of the pure 
Yang-Mills theory in ||3J), the hard thermal effective ac- 
tion for QCD, which can be obtained by applying the 
thermal operator, has a similar form to the one given in 
pTj) with the square of the thermal gluon mass given by 



g 2 T 2 



l QCD 



N - 



1 



-Nf 



(22) 



where Nf denotes the number of quark flavors. 



IV. CONCLUSION 

In this work, we have used the forward scattering de- 
scription for the retarded amplitudes in QCD to con- 
struct the hard effective action at one loop at zero tem- 
perature. By applying the thermal operator to the zero 
temperature amplitudes, we have derived in a simple way 
the hard thermal effective action for QCD. This approach 
emphasizes that various relevant features of this action 



5 



such as Lorentz invariance of the integrand in pTjl as well 
as its gauge invariance, simply arise because the leading 
hard zero temperature amplitudes (and the effective ac- 
tion) precisely have such properties. (The angular inte- 
gration breaks Lorentz invariance at finite temperature, 
as expected, since the rest frame of the heat bath defines 
a preferred reference frame. An alternative form with a 
Lorentz non-invariant integrand has been given in []|.) 
The above properties of the hard thermal effective action 
turn out to be very convenient in the study of the energy- 
momentum tensor for the quark-gluon plasma, as well as 
in the analysis of the high temperature behavior of gauge 
field theories in a curved space-time [lTJ. Furthermore, 
the effective action ([2T]) is also useful in implementing 
the resummation procedure [3j which is necessary for a 
consistent perturbative expansion in hard thermal QCD. 
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APPENDIX: HARD AMPLITUDES FROM WARD 
IDENTITIES 

For defmiteness, we discuss here in more detail some 
of the properties of the hard 3-point amplitudes at one 
loop. First, let us look at the ghost-ghost-gluon vertex, 
V^ hc (pi,P2, P3) at one loop and show how the leading 
terms cancel in this case. Let us examine the forward 
scattering graphs associated with such a vertex correction 
at one loop shown in Fig. 2, where p\ + p 2 + P3 = 0. 




FIG. 2: Examples of forward scattering graphs associated with the retarded ghost-ghost-gluon vertex at one loop. Color and 
Lorentz indices are suppressed for simplicity. 



In the hard regime (apart from an overall color factor 
f abc ) these graphs have a common numerator (p 2 ■ k)ka- 
The presence of quadratic terms in k in the numerator 
means that the leading terms, as in (2}, of the two re- 
tarded propagators (in the graph with k 2 = 0) are suf- 
ficient to give us the leading terms of degree zero in k. 
However, when the three graphs are summed, this con- 
tribution cancels as a consequence of the eikonal identity 



(k-p 1 )(k-p 2 ) (k-p 2 )(k-p 3 ) (k-p 3 )(k-pi) 

(A.I) 

Therefore, the ghost-ghost-gluon amplitude has a sub- 
leading contribution compared with the three point gluon 
amplitude, which as a result, satisfies the simpler Ward 
identity ©. 

Let us next show that the Ward identity (JSJ), together 
with the fact that to leading order the integrand of the 
3-point gluon amplitude, r ) c ^\(k,pi,p 2 ,p 3 ) 1 is a function 



of zero degree in k, is sufficient to determine the hard 
3-point gluon amplitude in terms of the hard gluon self- 
energy in We note that if we factor out an overall 
color factor of f abc , the Lorentz structure of the leading 
terms in (which are of degree zero in k and have the 
dimension of an inverse mass) can be parameterized in 
general as 

Bk^k v k\ + ^ {CxiPiykyky, + C 2t k^p^kx + C^k^Kpix) 

i 

+ (Eik^x + E 2 k u -q ll x + E 3 kxri^u) , (A. 2) 

where the coefficient functions B, Cu, C 2 i, C 3 i, E\, E 2 , E 3 
are Lorentz invariant functions of k and pi. The impor- 
tant point to note from (|A.2|) is that it is at most linear in 
Pi as well as in the metric tensor. (It is worth remarking 
here that while a Lorentz structure such as DijPi^pj U k\ 
is allowed in principle, for such a term to be of degree 
zero in k and have the inverse dimension of mass, the 
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coefficient Dij must contain a denominator of the form 
p 2 where p denotes a linear combination of the external 
momenta. However, as we have pointed out earlier, such 
terms do not arise in the expansion Q of the retarded 
propagator. This is why the numerator can at most be 
linear in the external momenta.) 

Given the general Lorentz structure (|A.2|) . let us next 
consider the Ward identity ([5]) by contracting the ampli- 
tude with p 3 which leads in the integrand to the Lorentz 
structure 

iPi -P3 k^k 

i 

+ ^2 {C 2l k ■ p 3 + EiSai) Piuk^ + E 3 k ■ p 3 j]^(A.3) 

i 

We can now compare this structure to the ones on the 
right hand side of the Ward identity in (jSJ coming from 



the structure of the self-energy in ([3]). Comparing the co- 
efficients of the metric tensor on both sides, the last term 
in (|A.3|) , namely, E 3 is uniquely determined. Likewise, 
the Ward identity © where we contract with or p^ 
determines respectively E\ and E2 uniquely. Similarly, 
matching the second and the third structures in (|A.3|) 
with the corresponding structures coming from the self- 
energy (along with the other two identities) uniquely de- 
termines J2iCiiPiv>I2i C 2iPiv'I2i C 3iPi\- Finally, com- 
paring the k^kv terms on both sides of the Ward identity 
([8]) then uniquely determines B. In this way, the leading 
3-point gluon amplitude is uniquely determined in terms 
of the hard gluon self-energy. One can follow this argu- 
ment to show that the leading higher point gluon am- 
plitudes are recursively determined uniquely by the hard 
gluon self-energy. This fact allows us to determine the 
effective action for gluon amplitudes in the hard region 
in a unique form. 



[1] J. I. Kapusta, Finite Temperature Field Theory (Cam- 
bridge University Press, Cambridge, England, 1989). 

[2] M. Le Bellac, Thermal Field Theory (Cambridge Univer- 
sity Press, Cambridge, England, 1996). 

[3] E. Braaten and R. Pisarski, Nucl. Phys. B337, 569 
(1990); ibidB339, 310 (1990); ibidPhys. Rev. D45, 1827 
(1992). 

[4] J. Frenkel and J. C. Taylor, Nucl. Phys. B334, 199 
(1990); ibidB374, 156 (1992); ibid B685, 393 (2004). 

[5] J. C. Taylor and S. M. Wong, Nucl. Phys. B346, 115 
(1990). 

[6] R. Jackiw and V. P. Nair, Phys. Rev. D48, 4991 (1993); 
R. Jackiw, Q. Liu and C. Lucchesi, Phys. Rev. D49, 6787 
(1994). 

[7] V. P. Nair, Phys. Rev. D48, 3432 (1993); ibid~D50, 4201 
(1994). 

[8] J. P. Blaizot and E. Iancu, Nucl. Phys. B390, 589 (1993); 
ibid B417, 608 (1994); ibid B434, 662 (1994); ibid Phys. 
Rep. 359, 355 (2002). 



[9] P. F. Kelly, Q. Liu, C. Lucchesi and C. Manuel, Phys. 
Rev. Lett. 72, 3461 (1994); ibid Phys. Rev. D50, 4209 
(1994). 

[10] A. Das and J. Frenkel, Phys. Rev. D75, 025021 (2007). 

[11] O. Espinosa, Phys. Rev. D71, 065009 (2005). 

[12] F. T. Brandt, A. Das, O. Espinosa, J. Frenkel and S. 

Perez, Phys. Rev. D72, 085006 (2005); ibidT>73, 065010 

(2006); ibid D73, 067702 (2006). 
[13] A. Das, Finite Temperature Field Theory (World Scien- 
tific, Singapore, 1997). 
[14] F. T. Brandt, A. Das, O. Espinosa, J. Frenkel and S. 

Perez, Phys. Rev. D74, 085006 (2006). 
[15] F. T. Brandt and J. Frenkel, Phys. Rev. D55, 7808 

(1997); ibidD56, 2453 (1997). 
[16] F. T. Brandt, A. Das, J. Frenkel and S. Perez, Phys. Rev. 

D74, 125005 (2006). 
[17] J. Frenkel and J. C. Taylor, Nucl. Phys. B410, 3 (1993); 

ibidB439, 131 (1995). 



